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ABSTRACT 

Context. Water fountain stars are very young post-AGB stars with high velocity water maser jets. They are the best objects to study 
the onset of bipolar jets from evolved stars due to their young dynamical ages. However, none of them have been observed in any 
thermal lines. 

Aims. Search for CO lines in the water fountain star IRAS 16342-3814 and investigate the properties of its thermal gas. 
Methods. The proximity, peculiar stellar velocity and high Galactic latitude of IRAS 16342-3814 make a single dish observation 
possible. We use the Arizona Radio Observatory 10m telescope to observe the CO J=2-l line and compare the line parameters with 
that of masers. 

Results. We report the detection of '-CO and "CO J=2-l lines from IRAS 16342-3814. The inferred '-CO mass loss rate is an order 
of magnitude lower than the infrared and OH mass loss rates, indicating a very cold and thick O-rich circumstellar envelope around 
the star. We also find a '^CO expansion velocity of Vexp =46+1 km s"' that is too high for an AGB wind and confirm the systemic 
velocity of 44 ± 1 km/s. In addition we measure a very low '^CO/'^CO line ratio of 1.7. 

Conclusions. The first detection of CO lines has provided a new way to investigate the water fountain stars. Given the high expansion 
velocity of the CO gas and its relation to maser velocities, we infer that the CO emission region is co-located with the OH mainline 
masers in the warm base of the optical bipolar lobes, while the high velocity OH 1612 MHz and H2O masers are located in the side 
walls and at the farthest ends of the bipolar lobes, respectively. Further observations are highly desired to understand the very low 
'^CO/'^CO line ratio. 

Key words, stars: AGB and post-AGB - (stars:) circumstellar matter - stars: evolution - stars: individual (IRAS 16342-3814) ~ stars: 
winds, outflows - radio lines: stars 

1 . Introduction which further constrained a higher line-of-sight H2O maser ex- 

pansion velocity of 130km s"'. 

Water fountain stars are a group of rare objects that display a Optical and infrared imaging to this star has revealed strik- 
thick circumstellar envelope (CSE) that is typical for asymp- ing structures in its CSE. T h e Hub ble Space Telescope (HST) 
C3 totic giant branch (AGB) stars and also high velocity bipolar jets observations by ISahai et alj (IT999h revealed a bipolar nebula. 
(> 100km s 1) detected via theii- H2O masers. They could be Corkscrew patterns were found in the bipolar lobes i n follow-up 
stars near the end of AGB or post-AGB evolution stage when infrared Lp band imaging with the Keck telescope jSahai et al.l 
non-spherical circumstellar structures begin to develop inside 2005), indicating that the bipolar lobes and the H2O masers 
the relic of the spherical AGB circumstellar envelope. Currendy might be created by p recessing jets. I nfrared imaging in the 
only about 10 such objects have been discovered through single 33^ 20yU range by Dii kstra et all ( l2003h showed the variation of 
dish and interferometry ob servations o f their OH, H2O and SiO the CSE morphologies from bipolar at near infrared to spherical 
masers (see the review by |lmaj|2007|). However, none of them at middle infrared, demonstrating that the outer CSE ejected ear- 
have ever been detected in a thermal molecular line. In this let- Her was spherical and th e star is current l y und ergoing CSE mor- 
ter, we report the first detection of the 7 = 2 - 1 transition of phology transformation. iDiikstra et all (|2003|) also found from 
i^CO and '^CO in the water fountain star IRA S 16342-3814. the 2 ~ 200// Infrared Space Observatory (ISO) spectra that the 
IR AS 16342-3814 was discovered by iLikkel & MorrisI star is the only object known to show absorption in both amor- 
(Il988h to show H2O masers with a line-of-sight expansion ve- phous and crystalline silicate features. The very red spectral en- 
locity up to 117.4kms"', much higher than the ty pical speed ergy distribution possibly indicates a very cold and thick dusty 
of~ 15kms ' for a spherical AGB wind (see, e.g.. lChen et alj CSE with an AGB mass loss rate on the order of lO^^Moyr ' 
l200lh . They also found a progressive increase of expansion ( IDiikstra et al.ll2003h . 

velocity in OH masers, from 36kms ' (1665 MHz) through Interferometric measurements of masers have supplied kine- 
42kms"' (1667 MHz) to 58 km s~' (1612 MHz). The symmet- matic inf ormation. The ma ppi ng of OH masers in I RAS 16342- 
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rical H2O maser velocities found bv lLikkel et an(ll992l) allowed 38 14 bv Sahai et alj (Il999t) and lZiilstra et alj ( I2OOII) showed lin- 



them to determine a star velocity of Vsys - 43.2 + 0.9 km s ear increase of line-of-sight outflow velocity with projected dis- 
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tance from the central star. High- velocity OH maser spots are lo- 
cated around the optical bipolar lobes while low-velocity maser 
spots appear near a dark waist close to the central star. These 
facts are in a ccord with a wind interaction model of the bipo- 
lar structure dZiilstra et al.ll200ll) . However, the corkscrew pat- 
tern in the infrared images also suggests that precessing jets 
may play some role. Very Long Baseline Array (VLBA) map- 
ping o f H2O masers by Morris et al. (20031) and IClaussen et al] 
( I2OO4I) showed bipolar distribution of maser spot clusters with 
more than 2" separation along the optical bipolar lobes and 
quasi-linear alignment of maser spots that delineates shock front 
of the jets. The interplay of the jets with magnetic field is 
still unclear, although extremely high degree of polarisation has 
been found in the OH masers dTe Lintel Hekkert et all [19881 
ISzvmczak & Gerardll2004 . 

The post- AGE n ature of IRAS 163 42-3814 has been criti- 
cally investigated by ISahai et al.l (11999 '). However the status of 
thermal gas in the CSE of this star is still unclear due to the 
lack of thermal molecular line observations. Unlike other known 
water fountain stars, IRAS 16342-3814 is the easiest object in 
which to search for its CO lines with a single dish telescope, 
because it is relatively near with a distance of 0.7 ~ 2 kpc 
( ISahai et al.lll999HZiilstra et al.ll2001h . it has a high galactic lat- 
itude of 5.°8 and a large line-of-sight velocity difference of at 
least 40 km s ' from that of nearby interstellar molecular clouds. 

2. Observation and data reduction 

The ALMA band-6 sideband separating receiver on the Arizona 
Radio Observatory 10 m submillimeter telescope (SMT) was 
used for the '^CO and '^COJ=2-l observations towards 
IRAS 16342-3814 on 2008, May 09. Sky subtraction was made 
in a beam switch mode with 2 arcmin throw at 2.2 Hz in azimuth 
direction. The receiver temperature was about 100 - 120 K. Two 
filter banks (FFBs, 1 GHz width, 1024 channels) were used for 
the two polarisations of the '^CO line and two acousto-optical 
spectrometers (AOS's, 972.8 MHz width, 2048 channels) were 
used for the two polarisations of the '^CO line. A 90 minute 
integration under a T225 = 0.15 and a T^ys = 445 K at a low 
elevation of 19° resulted in an RMS noise of about 10 mK at a 
resolution of 3 MHz in each polarisation. Main beam efficiencies 
of 0.7 and 0.6 were used to convert the measured into main 
beam temperature T^b for the USB ('^CO) and LSB ('^CO), re- 
spectively. A telescope efficiency of 3 1 Jy/K was used to derive 
line peak flux. The beam width is about 32" at both lines. 

The GILDAS/CLASS package was used to reduce and anal- 
yse the data. The two polarisations were combined to improve 
the S/N. A zero order baseline was removed from each spec- 
trum. We also experimented with a higher order baseline, but no 
much diff'erence was found at the band center where the CO lines 
appear. The SHELL method in the CLASS package was used to 
fit both lines with a truncated parabolic line profile. 

3. Results 

The resulting spectra of '^CO and '-^001=2-1 lines are shown 
in Fig.[Tl where each spectrum has been box smoothed to a res- 
olution of ~ 3 MHz. We do not see any strong interstellar emis- 
sion or absorption signatures. Line profile fitting of '^CO line 
gives a peak temperature of Tmb = 29 mK (with rms = 7 mK), 
corresponding to a peak flux of 0. 90 Jy (with rms = 0.22 Jy). The 
line area is 2.14 + 0.14Kkms"', representing a 15cr detection 
of the line. A systemic velocity of ysys('^CO) =44+1 kms"' 



and a high expansion velocity of yexp('^CO) = 46 + Ikms"' 
are measured by line profile fitting. The ysys('^CO) agrees well 
with the H2O mase r systemic velocity of 43.2 + 0.9 km s"' from 
iLikkel et al.l(ll992h . We also notice a pair of narrow features near 
the centre of the '^CO 2-1 line profile. However, because of the 
poor S/N, we are not sure if they are real features. 

The '^C02-l spectrum is more noisy. We find from line 
profile fitting a line peak temperature of rnib('^CO) = 17 mK 
(with rms - 9 mK), corresponding to a peak flux of 0.53 Jy (with 
rms = 0.28 Jy). The line area of 1.02 + 0.17 K km s"' represents 
a 6cr detection of the line. The stellar velocity and CSE expan- 
sion velocity found from line profile fitting are ysys('^CO) = 
50 + 3kms-' and V^^pC^CO) = 34 + 2kms"'. The velocity dif- 
ferences between the '^CO and '^C02-l lines could be merely 
caused by underestimated noise in the spectra. 

4. Discussion 

4.1. The low '^CO mass loss rate 

The detection of the '^CO 2- 1 line allows us to directly estimate 
the mass loss rate of IRAS 16 342-3814. Using the '^002-1 
formula scaled by Heske et al. (1990) from a '^CO 1-0 formula 
postulated by Knapp & Morris (1985.) , we derive M('-CO) = 
0.7 ~ 5.6 X 10~''Moyr-' for a source distance of 0.7 ~ 2 kpc. 
As discussed in § 14.21 the detected CO gas may not represent 
the spherical AGB wind. The true AGB wind CO mass loss rate, 
if also estimated with Heske's formula, should be even lower 
than this value, because the CO gas emission from the slow 
AGB wind is not definitely detected yet, whi le the mass loss 
rate is oc y^^p7'mb(C(9). However, the formula of lBaud & Habind 
( |1983h applied to an O H 1612MHz maser peak flux of 9Jy 
dLikkel&Morrisll 19881) yields a much higher mass loss rate of 
M(OH) = 1.7 ~ 4.9 X lO'^Mpy r'' (with V^„p = 46kms"' from 
our '^CO line). lDo et all (|2007|) fitted the Spitzer MIPS 70yU im- 
ages with a two shell model, and obtained M{lQfi) = 0.4 ~ 
3 X lO^^Moyr"' for an inner smooth CSE and M( shell) = 0.4 ~ 
3 X lO^^^Moyr"' for an outer extended dust shell (assuming a 
gas-to-dust mass ratio of 200 and Vg^p - 15kms"'. Here we 
note that our measured '^CO expansion velocity of 46kms"' 
may not be that of the dusty CSE, see the discussion in § 14. 2t . 
The inner shell mass loss rate agrees with the OH mass loss rate. 
Due to its huge radius of 1.4 ~ 4pc, the outer shell could be 
explained by sweep-up of interstellar mat erial, as su ggested by 
the recent hydrodynamic simulation of Wareing et al., (2007). To 
summarise, we suggest that IRAS 16342-3814, while located in 
a very extended and cold detached outer dust shell, has an inner 
CSE whose mass loss rate (estimated from IR and OH emission) 
is roughly an order of magnitude higher than its '^CO mass loss 
rate. The weakness of '^CO 2-1 line could be interpreted by the 
coldness of th e CO gas, as sugge sted for some other very red 
OH/IR stars bv lHeske et akl d 19901) . 

4.2. The high '^CO expansion velocity 

The '^CO expansion velocity yexp('^CO) - 46kms"' is too 
high for a typical OH/IR star (typ ically ~ 15 + 5 kms~' for 
long period OH/IR stars, see e.g.. IChen et al.ll2001i) . Another 
possible source of the broad '^CO emission is the swept-up 
bipolar bubble walls. According to a momentum driven wind- 
interaction model suggested by Ziilstra et al. (2001), the radial 
motion velocity of different parts of the bubble wall will be pro- 
portional to distance to the central star This has been confirmed 
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by the li near OH maser vel ocit y ofFset-radius (V-r) relationship 
found bv lSahai et all d 19991) and lZiilstra eTall (1200 lb . Assuming 
that VexpC'^CO) fo llows the same V-r relationship in Fig. 3 of 
ISahai et al.l (199^ . we find a rough angular radius of 0.5" for 
the CO emission region, similar to that of OH mainline masers. 
Further assuming a Gaussian source shape, a brightness tem- 
perature of the '^C02-l line can be estimated to be ~ 30 K. 
Therefore, the '^C02-l line could be optically thin if the gas 
temperature is much higher than this. 

Further insight can be gained by comparing the velocities 
of '^C02-l line, OH and H2O masers. We find several inter- 
esting facts by comparing the CO velocit ies with the V-r re - 
lationship of OH masers in the Fig. 3 of ISahai et a D (Il999h : 
1) The 1612MHz maser velocity offsets are all > VexpC^^CO), 
while most of the OH 1665/1667 MHz (mainline) maser veloc- 
ity offsets < VexpC'^CO); 2) The OH 1612 MHz masers show 
two components - a pair of bipolar high velocity clumps that 
tightly follow the linear V-r relation and a peculiar group of 
points with Vlsr ~ Okms"' but scattered in a large radius 
range (-1.1" ~ 0"); 3) The position of the bipolar 1612MHz 
maser component is not symmetrical with respect to the star po- 
sition, perhaps due to 3D projection effect. The peculiar group of 
1612 MHz masers are not only scattered in radius, but also scat- 
tered in position angle with resp ect to the central star in the VLA 
maps (see Fig. lb and Fig. 9 of Sahai et alj|l999t IZiilstra et alj 
[2001, respectively). This is difficult to explain with a spherical 
expanding shell model in which masers should show an ellipti- 
cal relation between line-of-sight velocity and projected radial 
distance. Taking away this peculiar 1612 MHz maser compo- 
nent, one can see a better V-r relationship among the rest 1612, 
1665 and 1667 MHz masers. If we further take a velocity span 
of ~ 248kms~' and a s patial separation of ~ 3" of H2O masers 
from lMorris et al.l (12003 ) and overlap their symmetry center with 
that of OH masers, we can see that the H2O masers follow the 
same linear V-r relationship quite well. Upon these arguments, 
we can sketch a picture for the CSE as follows: 

Both the OH mainline masers and the CO emission arise 
from the same region near the base of the bipolar outflow, the 
two clumps of high velocity OH 1612 maser spots are located 
on the farther and faster part of the bipolar lobe side walls, the 
H2O masers are stimulated by the shock fronts at the two farthest 
and fastest ends of the bipolar lobes, while the peculiar group of 
spatially scattered 1612 MHz maser spots are located in an in- 
dependent structure that is blue-shifted by ~ -46kms"' with 
respect to the star but not necessarily associated with the bipolar 
outflow. The dark waist and the slow AGB wind might be simply 
too cold to be detected in our CO 2-1 line spectra. Actually, the 
pair of faint spikes on the top of the broad '^CO 2-1 line profile 
near the line centre in Fig.[T]could be the emission from the slow 
AGB wind. If this is true, we may estimate a crude AGB wind 
speed of ~ 13 km s \ typical for an AGB star 

4.3. On the origin of the extremely low ^^COA^CO ratio 

The detection of the '^CO J=2-l line yields an extremely 
smafl '^CO/'^CO peak flux ratio of 1.7. A range of '^CO/'^CO 
ratios ha ve been fo und in PPNe and PNe by other re- 
searchers. iHasegawal (l2005h con cluded frorn a siirvey of pub- 
lished CO observations o f PNe dPall a et al."2000': 'Bal ser et alJ 
I2OO2I: lBachilleretaLlll997l: iJosselin & Bachifler 2003) that the 
^^CO/^^CO ratios are about 10 ~ 30 in most PNe, but are 
2 ~ 3 in a few and > 60 in a few others. Some ex- 
amples of low ratios are: ^^CO/ ' ^CO ~ 4.6 in AFGL618, 
~ 2.2 in M 1-16 dBalser et al.l l2002h . ~ 3 in the PN 



NGC6302 (pasegawa & Kwoai2003h. and - 3.2 in HD 17982 
(iJossehn & Lebrell2001l ). 'Bachill er etall d 19971) pointed out that 
the '^CO/'^CO ratios in PPNe and PNe are affected by two com- 
peting effects: selective photodissociation of '^CO and isotopic 
fractionation of '^C in CO. Another issue is the opacity in the 
'^CO lines, e.g., the very l ow ratio in NGC 630 2 was later found 
by interferometry work of Perett o et al.l d2007h to be affected by 
opacity effects in a massive torus and the true '^C/'^C should be 
> 15. For an embedded dusty post- AGB star like IRAS 16342- 
3814, we do not need to worry about the selective photodisso- 
ciation in the central region. If the detected CO emission does 
come from the accelerated bipolar bubble walls (see the discus- 
sion in § |4.2t , neither isotopic fractionation nor '-CO line opac- 
ity should be an issue, since the dynamical age of the bipolar 
jets is short (~ 100 yr) and the high expansion velocity and high 
radial velocity gradient (according to wind interaction model) 
might render the '^CO emission not too optically thick. A rough 
but conservative calculation of the optical depth under an LVG 
approximation gives a value smaller than unity. However, better 
CO line data is required to confirm this. We also find that our 
'^CO/'^CO ratio is about 1.1 and 3.2 on the red and blue line 
wings (46 > IVexpl > 15kms"') respectively, but the uncertain- 
ties are larger there. One more issue is that the '-CO and '^CO 
lines might come from different clumps of gas. Again, better CO 
line profiles and, better still, CO mapping data are needed to 
clarify this point and thus to constrain a reliable '^C/'^C ratio. 
Therefore, later in this section, we will discuss different opacity 
cases separately, using only the assumption that the '^CO and 
'^CO lines are co-located in the CSE. 

A precise measurement of the '^C/'^^C ratio has the poten- 
tial to cast light on the evolutionary status of the star. Stellar 
models of low mass (M < 2.0Mq) and intermediat e mass 
(2.OM0 < M < 8.OM0) stars (e.j 



, |E1 Fid 1994; Stranier o et al.l 
ll997l:lKarakas & Lattanziol2007l) show that there are quite a few 
phases of evolution where this ratio is altered from the initial 
value (assumed to be solar at 89). The first occurs during the 
first dredge up (FDU) at the beginning of the Red Giant Branch 
(RGB). Here the ratio declines from the initial value to ~ 20 as 
the convective envelope moves in and mixes up partially burnt 
material (we note that the second dredge-up has only a minor 
effect on the ratio). The second evolutionary phase that alters 
the '^C/'^C ratio is not predicted by stellar models. Observations 
have shown that the isotopic ratio reduces along the RGB in low- 
mass stars (M < 2.OM0), via a process known as Deep Mixing 
(als o known as 6u m ixing - see the recent theoretical models by 
Egg leton et al.ll2008l) . This phenomenon reduces the ratio down 
to a value of ~ 15 by the end of the RGB in Population I stars. 
The third phase that significantly alters the ratio actually in- 
creases it - the Third Dredge-Up (TDU). Here '^C is periodi- 
cally mixed into the envelope during the thermally-pulsing AGB 
stage (TPAGB), raising the '^C/'^C r atio from the FPU v alue 
to values of ~ 60 or more (see e.g.. iLebzelter et ani2008l for 
recent observations of this effect). The final phase to alter the ra- 
tio is Hot Bottom Burning (HBB). However this only occurs in 
stars more massive than ~ 4 or 5M0 (at solar metallicity). HBB 
is active during the later phases of the TPAGB and acts to burn 
the carbon and nitrogen in the convective envelope via the CNO 
cycles. This results in the '^ C/'-^C ra tio appro aching the equi- 
libriu m value of ~ 3 (see e.g.. lPlez et a l. 1993; McSavenev et al.l 
l2007',for some observational evidence of this phenomenon). The 
very low '^C/'^C ratio in some Galactic PNe, PPNs and ex- 
tremely young post- AGB stars like IRAS 16342-3814 could be 
also indicative of the HBB process. 
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If IRAS 16342-3814 has a high '^C/'^C ratio close to the 
solar value of 89 (i.e., ~ 50 times larger than the '^CO/'^CO line 
ratio, e.g., as a 3 star experiencing TDU), the CO line optical 
depths can be derived from the observed '^CO/'^CO ratio (=1.7) 
to be t(12C0) ~ 82 and t('^CO) ~ 0.93 (according to T oc 
(1 - e"^) under the assumption of co-location of '^CO and '^'CO 
gas and the same excitation temperature). However, the '^CO 
mass loss rate should also be raised by a factor of r/( 1 -e^'^) = 82 
to 5.7 ~ 45.9 X lO^^Mo/yr, which is significantly higher than the 
OH and IR mass loss rates. 

If the true '^C/'^C ratio is an intermediate value of 20, i.e., a 
factor of ~ 10 higher than the '^CO/'^CO line ratio, the similarly 
determined opacities are t(^^CO) ~ 19 and t('^CO) ~ 0.93. In 
this case, the '^CO mass loss rate is comparable to the OH and IR 
mass loss rates. Therefore, a low mass AGB star (< 1.5 Mq) that 
didn't experience TDU or HBB and thus retains its FDU plus 
Deep Mixing surface abundances would be a possible progenitor 
of IRAS 16342-3814. 

If the '^C02-l line is only slightly optically thick so that 
the '^C/'^C ratio is ~ 3 (the CNO equilibrium value), the CO 
Une opacities are t^^CO) ~ 2.3 and r('^CO) ~ 0.77. In this 
case, we can see that the only process (to the best of our knowl- 
edge) that could have lowered the '^C/'^C ratio to such a value 
is HBB. Thus the inferred mass of the star that produced the 
material that we observe today would have been > 4 or 5 Mq. 
Interestingly the low '^CO/'^CO ratio, bipolar outflow, and the 
properties of the OH masers of IRAS 16342-3814 are remi- 
niscent of the ma s sive P PN IRAS 220364-5306 investigated by 
ISahaiet al.1 (I2003L l2006h . although IRAS 16342-3814 is proba- 
bly in an earlier phase of evolution due to the presence of the 
young 'water fountain'. Alternatively, a binary system in which 
a massive AGB star had transfered low '^C/'^C matter to the 
current AGB star is also possible. We note that a low isotopic 
ratio is supported by our rough estimation of the small '^CO 2-1 
opacity at the beginning of this section. 

However, both the measured CO and H2O maser systemic 
velocities (~44km s '. lLikkel et al.|[T992l) confirm that the spa- 
tial motion of IRAS 16342-38 14 is in the reverse direction of the 
Galactic disk rotation (interstellar CO cloud ve locity near the 
same line of sight is ~ -17kms"', Dame et al. 2001). The pe- 
culiar stellar motion lends weight to IRAS 16342-3814 being a 
low mass stars, although in rare cases kick-out of stars from the 
crowded galactic plane can give rise to peculiar velocities. 

Given the preceeding discussions it is obvious that further 
observations are highly desirable to investigate the different 
opacity possibilities. Molecular line interferometry and maser 
prope r motion measurements (such as the work of llmai et alj 
12007 ) are also needed to check the interstellar contamination, 
confirm the co-location of the '^CO and '^CO gas, and to locate 
the birthplace of the possibly escaped AGB progenitor on the 
Galactic plane. 
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. 1. The observed '^CO and '^CO 2-1 line profiles. The spectra have been box smoothed to a resolution of about 3 MHz. The 
O 2- 1 line has been shifted upward by 0.06 K for clarity. The arrows mark Vsys = 44 km s" ' that is derived from the '^CO Une. 



